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ABSTRACT 



We report the discovery of a "super-damped" Lyman-alpha absorber at Zabs — 
2.2068 toward QSO Ql 135-0010 in the Sloan Digital Sky Survey, and follow- 
up VLT UVES spectroscopy. Voigt profile fit to the DLA hue indicates log 
A^Hi = 22.05±0.1. This is the second QSO DLA discovered to date with such high 
Ahi. We derive element abundances [Si/H] = -l.lOiO.lO, [Zn/H] = -1.06±0.10, 
[Cr/H] = -1.55 ±0.10, [Ni/H] = -1.60 ±0.10, [Fe/H] = -1.76 ±0.10, [Ti/H] = 
-1.69±0.11, [P/H] = -0.93±0.23, and [Cu/H] = -0.75±0.14. Our data indicate 
detection of Ly-a emission in the DLA trough, implying a star formation rate of 
~10 Mq yr~^ in the absence of dust attenuation. C II* A1336 absorption is also 
detected, suggesting SFR surface density —2 < log^'* < Mq yr'^ kpc~^. We 
estimate electron density in the range 3.5 x 10~^ to 24.7 cm~^ from C II*/C II, 
and ~0.5-0.9 cm~^ from Si II*/Si II. Overall, this is a robustly star-forming, 
moderately enriched absorber, but with relatively low dust depletion. Fitting of 
the SDSS spectrum yields low reddening for Milky Way, LMC, or SMC extinction 
curves. No CO absorption is detected, and C I absorption is weak. The low dust 
and molecular content, reminiscent of some SMC sight-lines, may result from the 
lower metallicity, and a stronger radiation field (due to higher SFR). Finally, we 
compare this absorber with other QSO and GRB DLAs. 

Subject headings: Quasars: absorption lines-ISM: dust 
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1. Introduction 

A flux-independent probe of galaxy evolution is provided by absorption lines in spectra 
of quasi-stellar objects (QSOs) superposed by foreground galaxies along the sightlines, 
particularly the damped Lyman-a (DLA; neutral hydrogen column densities A^hi > 2 x 10^° 
cm"^) and sub-DLA (10^^ < Ahi < 2 x 10^° cm"^) absorbers. DLAs and sub-DLAs are the 
primary neutral gas reservoir available for star formation (e.g., Storrie-Lombardi & Wolfe 
2000; Peroux et al. 2005; Prochaska et al. 2005). Furthermore, they offer among the most 
accurate element abundance measurements in distant galaxies. The neutral hydrogen and 
metal content of DLAs/sub-DLAs provide unique indicators of the star formation history, 
regardless of galaxy redshift or morphology. 

The calculation of the co-moving density of neutral hydrogen {Q^hi) directly relies on 
the integration of the QSO absorber Ahi column density distribution, f{Niii). The slope 
of /(Ahi) at large Ahi is typically a = —3.5, implying that systems with large column 
density are very rare. Noterdaeme et al. (2009) reported a definite detection from UVES 
high spectral resolution data of the first DLA with log A^hi = 22.0 ± 0.1, at Zabs — 3.286 
towards SDSS J081634+144612. This was the highest H I column density DLA reported so 
far along a QSO line of sight. Such a large H I column density is similar to that of DLAs 
detected at the redshifts of some Gamma-Ray Bursts (GRBs; e.g. Watson et al. 2006; 
Ledoux et al., 2009). Existence of such strong QSO DLAs puts a strong constraint on the 
shape of the H I column density distribution. Indeed, from the fit of the H I column density 
frequency distribution, Noterdaeme et al. (2009) argue that no more than one system like 
this is expected in the whole SDSS survey. 

^Based on observations collected during program ESO 385.A-0778 at the European South- 
ern Observatory with UVES on the 8.2 m KUEYUN telescope operated at the Paranal 
Observatory, Chile. 



-4 - 



Recently, we serendipitously discovered another such high A^hi DLA at Zabs — 2.2068 in 
the spectrum of QSO SDSS J113520.39-001053.5 with z^ra = 2.907 (see Fig. 1). The QSO 
is not a broad-absorption- hne QSO. The SDSS spectrum indicates that the DLA has log 
A^Hi ~ 21.9. This system is probably missing from the sample of Noterdaeme et al. (2009) 
because it is outside their automated statistical search which has a fixed minimum redshift 
per spectrum corresponding to a threshold in S/N. The damped Lyman-alpha line is in 
fact not completely covered in the SDSS, and the system falls below the low- redshift edge 
of the redshift sensitivity function of both Noterdaeme et al. (2009) and Prochaska & Wolfe 
(2009). 

Several strong heavy element lines are seen in the SDSS spectrum which suggest a 
relatively metal-rich absorber. This DLA is thus an ideal place to search for rare elements 
and molecules. Of course, to obtain accurate abundances, higher resolution spectra are 
essential. Here we report follow-up VLT UVES spectroscopy of this unique sighthne. 

2. Observations and Data Reduction 

SDSS J113520.39-001053.5 was observed in queue mode on April 13, 2010 with VLT-U2 
using UVES (Dekker et al. 2000). Observing time of about 16,200 s was awarded for 
this program, but integrations of only 4,760 s were actually obtained due to scheduling 
constraints. The average seeing was 0.99 Two settings centered at Ac = 390 nm and 
Ac = 564 nm were used to cover observed wavelengths in the range of about 3290-6648 
A (with some gaps in the ranges of 4518-4625, 5600-5674). This corresponds to a DLA 
rest-frame coverage of 1026-1409, 1442-1746 and 1769- 2073A, giving access to key lines of 
many metal ions. A 1" slit was used, with 2x2 binning. 

The data were analyzed with the UVES pipeline, IDL and IRAF. The red side was 
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re-processed with the latest reflex pipeline. The blue side was binned down to ~ 6000 
assuming a 3 pixel resolution element, to improve the S/N. This resolution is adequate for 
fitting the damped Lyman-alpha line. The red side data were binned by a factor of 2, as 
these spectra, taken with a 1" sht and 1x1 binning, were well oversampled. The binned 
spectra have a resolution R ~ 45000 per 3 pixel resolution element. The continuum levels 
were determined using cubic spline fits to ~ 20 A-wide individual sections of spectra. The 
absorption lines were fitted with Voigt profiles using the program FITS6P that we have used 
extensively in previous studies of interstellar and intergalactic matter. FITS6P (Welty et 
al. 1991) evolved from the code used by Vidal-Madjar et al. (1977), and minimizes the 
between the data and the theoretical Voigt profiles convolved with the instrumental profile. 



3. Results 

Fig. 2 shows the damped Lyman-alpha absorption profile covered in our UVES blue 
spectrum. Fig. 3 and 4 show velocity plots of key absorption lines of various ions. Voigt 
profile fitting gives H I column density of log Nui = 22.05 ± 0.10. The solid red and dashed 
blue curves overlaid on the data in the bottom panel of Fig. 2 show the Voigt profile fits 
corresponding to log A^hi = 22.05 ± 0.10. The upward points in the region 4030-4040 A are 
near the location of the Ly-/3 and O VI emission hues of the QSO, and also comparable 
to the noise level (the latter is shown in the top panel of Fig. 2 as a dashed black curve). 
It is clear that the absorber has a very high A^hi- For comparison, we have also shown in 
dot-dashed orange curves the Voigt profiles for log A^hi = 21.85 and 22.25, which are clearly 
too extreme. Indeed, log A^hi = 21.85 is too low to explain the strong core of the line even in 
the region ~ 3920 - 3960 A which appears relatively free from blends with other lines. We 
adopt log A^Hi = 22.05 ± 0.1 as our best estimate of the H I column density. We use this A^hi 
estimate in deriving the element abundances below, since the contribution of ionized gas is 
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not expected to be large at such a high H I column density, and since the H2 contribution 
is hard to determine (see section 4.1.1). We note that the DLA is not a superposition of 
two smaller column density clouds, which would otherwise have been reflected in the metal 
lines as well. 

3.1. Absolute Abundances 

The spectra show lines of many ions-low ions such as C I, C II, C IV, O I, Al II, Al 
III, Si II, Si IV, P II, Ti II, Cr II, Fe II, Ni II, Cu II, and Zn II. The multi-component 
Voigt proflle flts to the metal lines are shown overlaid on the data in Fig. 3 and 4. Tables 
1 and 2 list the column densities derived for individual components. For each component, 
we have hsted the mean velocity relative to 2; = 2.2068, the Doppler b parameter, and the 
column density. The element abundances summed over all the velocity components were 
calculated with reference to solar levels from Lodders (2003), and are listed in Table 3. 
We infer [Zn/H] = -1.06 ± 0.10, [Fe/H] = -1.76 ±0.10, [Cr/H] = -1.55 ± 0.10, [Ni/H] = 
-1.60 ±0.10, [Cu/H] = -0.75 ±0.14, [P/H] = -0.93 ± 0.23, [Ti/H] = -1.69 ± 0.11 and 
[Si/H] = —1.10 ± 0.10. Overall, this is a moderately enriched absorber. Some of the column 
density measurements have substantial uncertainties, especially those for P and Cu. We 
note, however, that, to the best of our knowledge, this is the first detection of Cu in a DLA. 

3.2. Relative Abundances 

To examine the nature of dust depletion, we calculated the relative abundances [Fe/Zn], 
[Cr/Zn], [Ni/Zn], [Cu/Zn], [Cu/Fe], [Fe/Si], [Ti/Si], and [P/Si] in the individual velocity 
components, and compared them to the interstellar warm cloud depletion values for the 
Milky Way, based on a combination of Welty et al. (1999), Cartledge et al. (2006), and 
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Jenkins (2009). We find that there is considerable variation in the relative abundances from 
component to component, indicating variation in dust content. Surprisingly, Ti is the most 
depleted element in only two of the seven velocity components that show Ti detections, 
but is less depleted than Fe, Ni, and/or Cr in the remaining components, with Ni and Cr 
showing smaller abundances than Fe relative to Zn by 0.2-0.5 dex in all but one component. 
It is also interesting to note that the abundance of Ni relative to Cr is inverted in the 
component at v = 62.5 km s~^. The relative differences between Ti and Fe, Cr, Ni are 
reminiscent of some velocity components in the SMC interstellar sighthne toward Sk 155 
(e.g. the [Ti/Zn] variations in Welty & Crowther 2010). For most abundance ratios such as 
[Fe/Zn], [Cr/Zn], [Ni/Zn], [Fe/Si], the values observed in the DLA under study are higher 
(indicating less severe depletions) than in the Milky Way warm diffuse clouds, by differing 
amounts in different velocity components. However, [Cu/Fe], [Ti/Si], and [P/Si] are below 
the warm diffuse cloud levels for some velocity components. This may suggest different dust 
grain compositions in the super-DLA. On the other hand, the low [Cu/Fe] and [P/Si] values 
in some components may be intrinsic nucleosynthetic differences, reflecting the odd-even 
effect. 

Fig. 5 shows the relative abundances as a function of the Zn II column density in 
the velocity components with detectable Zn. The Pearson product-moment correlation 
coefficient for [Ni/Zn] vs. log A^znii is -0.961, which suggests a fairly strong anti-correlation 
(the 1-tailed probabihty of the correlation happening by pure chance being 0.0003). For 
[Fe/Zn] vs. log A^znii, the Pearson product-moment correlation coefficient is -0.455, which 
implies no significant correlation (1-tailed probability = 0.152). For [Ti/Si], the coefficient 
is -0.823, with a 1-tailed probability of 0.0221 of this correlation arising purely by chance. 
Thus, the relative abundances of some elements are consistent with being correlated with 
the amount of metals (traced by the Zn II column density). While this trend is far from 
robust (since it pertains to dust in a small number of velocity components of one absorber). 
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it is reminiscent of the total [Cr/Zn] vs. total metallicity correlation reported for DLAs 
and sub-DLAs (e.g., Ledoux et al. 2003, Meiring et al. 2009a,b and references therein). We 
note, however, that such correlations cannot always be explained by dust, because similar 
correlations are also observed when the relative abundance does not indicate dust depletion, 
e.g., for [S/Zn] vs log A^znii in both QSO DLAs and GRB DLAs (e.g., Savagho et al. 2012). 

Fig. 6 shows the total abundances (summed over all velocity components) of Cu, P, Si, 
Fe, Ni, and Ti relative to Zn, plotted as a function of the condensation temperature of the 
element (the latter adopted from Lodders 2003). For each element, the black circles denote 
the data points for the DLA toward Q1135-0100. The red squares and blue diamonds 
denote the relative abundances for warm and cold diffuse interstellar clouds in the Milky 
Way. Elements with higher condensation temperatures normally show higher depletions in 
the Milky Way clouds, with warm clouds showing smaller degree of depletions than cold 
clouds. The relative abundances [P/Zn] and [Si/Zn] in the DLA toward Q1135-0100 are 
close to zero, at levels expected for warm Milky Way clouds. The relative abundances 
[Cr/Zn], [Fe/Zn], [Ni/Zn], and [Ti/Zn] are higher, as expected for the Milky Way clouds, 
but are less severe than even those in the warm clouds. For Cu, there is a large excess, with 
[Cu/Zn] being super-solar at about the 2 a level. This is surprising, since both depletion 
and intrinsic nucleosynthesis (odd-even effect) would predict [Cu/Zn] < 0. 

3.3. Other Ions 

One signature of a very gas-rich DLA is the presence of weak C I lines. In the DLA 
toward Q1135-0010, CIA 1656 is detected, although in only the velocity component at -39 
km s^^. Voigt profile fitting gives A^(CI) = (8.89 ± 1.51) x 10^^ ^he latter is consistent 
with the non-detection of C I A 1560. The observed ratio of C 1/ H I is low for the high H I 
column density, compared to the Milky Way, and may be a result of the lower metallicity 
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and a stronger radiation field. 

A few liiglier ionization species are also detected in this absorption system. Table 4 
lists the column densities for these estimated using Voigt profile fitting. Fig. 7 shows a 
velocity plot of these ions. 

3.4. Lyman-o; Emission 

Despite the noise in our spectrum, there is a hint of Lj-a emission near the center of the 
DLA absorption trough. To better assess the strength of this feature, we resampled the data 
by varying factors and found the feature present in all cases. Fig. 8a shows the UVES data 
resampled by factors of 30 and 70, along with the SDSS data. Fig. 8b shows the UVES data 
resampled by a factor of about 70 (to a resolution of 1 A) , superposed with the same data 
box-car smoothed by factors of 5, 11, and 21. It is clear that the feature does not go away. 
In the original SDSS spectrum, the integrated flux in the feature is (2.87 ± 0.63) x 10~^^ erg 
s^^ cm^^, where the flux uncertainty is estimated as the average photon noise uncertainty 
outside the feature at the bottom of the DLA trough, multiplied by the square root of the 
number of pixels covered by the feature. In the UVES spectrum resampled by factors of 
30, 50, and 70 (i.e., to resolutions of 0.42, 0.70, and 0.98 A respectively), the integrated 
counts in the feature are 29.1 ± 10.9, 46.3 ± 12.2, and 77.0 ± 16.0, respectively, indicating 
detections at levels of 2.7, 3.8, and 4.6(7 levels. [The resultant estimates of the integrated 
flux in the Ly-a emission line would be (1.15 ± 0.45) x 10"^^ (1.15 ± 0.30) x 10"^^ 
and (1.11 ± 0.20) x 10^^^ erg s^^ cm^^ respectively, if the flux calibration of the UVES 
spectrum were used. The difference between the SDSS and VLT line flux values arises due 
to differences in the detailed structure of the observed Ly-a emission features in the SDSS 
and VLT data, and differences in the absolute flux cahbration of the SDSS and VLT data.] 
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We use the SDSS spectrum to estimate the absolute flux in the Ly-a emission feature, 
since the SDSS absolute photometric calibration is believed to be much more accurate. 
Thus we adopt the integrated Ly-a emission flux to be (2.87 ± 0.63) x 10~^^ erg s~^ cm~^. 
Assuming Hq = 70 km s~^ Mpc~^, Qm — 0.3, and Qa — 0.7, the luminosity distance at 
z = 2.2068 is 17.522 Gpc. We thus estimate a Ly-a luminosity of 1.05 x lO''^ erg s'^, 
i.e. 2.7 X IO^Lq. Assuming case-B recombination and adopting the Kennicutt (1998) 
relation between the star formation rate (SFR) and R-a luminosity, we estimate SFR = 
9.1 X 10~'^^ Liy-a — 9.6 ± 2.1 Mq yr~^. These estimates assume no correction for dust 
absorption or resonant scattering by neutral gas. As we discuss in section 4.3, there is not 
much evidence of a signiflcant amount of dust in this DLA. The inferred Ly-a emission and 
SFR, if conflrmed with higher S/N spectra in the future, is much stronger than in most 
other QSO absorbers (e.g., Bouche et al. 2007; Fynbo et al. 2011; Kulkarni et al. 2000, 
2006; Rahmani et al. 2010; Peroux et al. 2011a, 2011b, and references therein). Indeed, as 
noted by Wolfe & Chen (2006), the star- formation efficiency for high-redshift DLAs is low. 
We discuss the star formation in the DLA toward Q1135-0100 further in section 4.4 using 
C II* absorption. 



4. Discussion 

4.1. Comparison with Other High-iVni Absorbers 

We now compare the super-DLA toward Ql 135-0100 with other DLAs observed toward 
QSOs and GRBs. The comparison QSO DLA sample is from the compilation used in 
Kulkarni et al. (2010), while the GRB DLA sample is based on Prochaska et al. (2007, 
2009), Ledoux et al. (2009), Rau et al. (2010), and Savaglio et al. (2012). For ease of 
comparison, we have only used systems with firm measurements of Zn (in most cases) or S 
(in a few cases with no Zn detections), excluding upper or lower limits. Both Zn and S are 
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almost undepleted in the warm Milky Way interstellar clouds, and can thus serve as good 
indicators of the total metallicity, without the need to correct for dust depletion. 

4.1.1. Metallicity , H I Content, and Molecular Gas 

The possibility that some high-redshift QSOs may be obscured by dust in foreground 
DLAs was suggested by Fall & Pei (1993). Boisse et al. (1998) pointed out that there is a 
deficit of DLAs with high A^hi and high metaUicity. Since the metal lines of such systems 
should be easy to measure, Boisse et al. postulated that the deficit is a selection effect 
caused by dust obscuration: systems with large H I column densities and high metallicities 
would have larger metal column densities and hence larger dust column densities. Such 
absorbers could therefore obscure the background QSOs enough that they would tend to 
be underrepresented in spectroscopic studies that often tend to favor a brighter QSO over 
a fainter QSO. Boisse et al. suggested an empirical obscuration threshold of log A^znii = 
13.15. Over the past few years, a few DLAs above this threshold have been discovered 
toward QSOs and GRBs (e.g., Meiring et al. 2006; Peroux et al. 2008), thereby casting 
doubt on the empirically inferred existence of the obscuration threshold. Several of the 
systems detected beyond the threshold, however, lie close to the threshold. 

Fig. 9 shows a plot of the metallicity vs. H I column density for QSO DLAs and 
GRB DLAs. The QSO DLAs are from Kulkarni et al. (2010) and references therein. 
The GRB DLAs are from Vreeswijk et al. (2004), Prochaska et al. (2007), Ledoux et al. 
(2009), Savaglio et al. (2012), and references therein. The solid orange line denotes the 
obscuration threshold suggested by Boisse et al. (1998). The DLA toward Q1135-0100 is 
significantly beyond the "obscuration threshold"; in fact, it is the "farthest" QSO DLA from 
the obscuration threshold known so far. In terms of the "distance" from the obscuration 
threshold line, the DLA toward Q1135-0100 is in fact, comparable, to the DLA toward 
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GRB 000926. The DLA toward GRB080607 is believed to be even further beyond in the 
"zone of avoidance", although its metallicity is highly uncertain owing to line saturation 
(Prochaska et al. 2009). 

A possible explanation of the low A^hi in QSO-DLAs was suggested by Pontzen et al. 
(2010): QSO-DLAs are in random sight lines, likely in the external part of galaxies where 
most of the volume is and where gas density and metaUicity are lower, whereas GRB-DLAs 
are in star-forming regions closes to the galaxy center, where gas column densities are on 
average larger. However, this does not explain why the metallicities of most QSO DLAs 
are lower than those in QSO sub-DLAs (e.g., Kulkarni et al. 2010 and references therein), 
which should be tracing gas even further out. 

Another explanation of the deficit of DLAs with high Ahi and high metallicity was 
suggested by Schaye et al. (2001) and Krumholz et al. (2009a), in terms of the formation of 
molecules in the cold phase of a 2-phase interstellar medium in pressure balance. Their idea 
is that above some Ahi threshold that decreases with increasing metallicity, the hydrogen 
becomes predominantly molecular, and hence undetectable in observations of the Ly-a 
absorption line. A fixed optical depth of material can absorb the entire photon flux in the 
Lyman/ Werner region, so any additional gas is molecular. A higher metallicity implies a 
higher dust-to-gas ratio, which results in an increased H2 formation rate. This, because of 
the increased cooling via metals and H2, result in a lower temperature. Both the higher H2 
formation rate and the lower temperature lead to a decrease in the maximum A^hi- 

The remaining lines in Fig. 9 show curves calculated by Krumholz et al. (2009a) for 
"covering" fractions c^^ = 0.01, 0.05, 0.5, and 1.0 for the cross-section of the spherical 
molecular core of the cloud (surrounded by a shell of atomic gas). These values of 
correspond to molecular mass fractions /hj = -^HaZ-^totai — / (l^moi}]'^ — 0.010, 

0.102, 0.845, and 1.0, respectively. [The quantity (f)moi is the ratio of molecular gas density 
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to atomic gas density, and is about 10 (Krumholz, McKee, & Tumlinson 2009b)]. In this 
simple picture, the DLA toward Q1135-0100 may seem to be consistent with an H2 core 
"covering" fraction of > 0.05 (although values as small as ~ 0.03 and as high as ~ 0.5 
cannot be ruled out within the uncertainties on the metallicity and H I column density). 
Thus the core "covering" fraction could be higher than the values inferred to be < 0.06 
by Krumholz et al. (2009a) for most other QSO or GRB DLAs (except that toward 
GRB080607; Prochaska et al. 2009). 

In this scenario, the molecular mass fraction for the DLA toward Q 1135-00 10 could 
be as low as 0.05 or as high as 0.85-much higher than that observed for other QSO DLAs, 
in either case. This DLA would therefore be a good candidate to search for molecules. 
Unfortunately, the Lyman and Werner bands of H2 in the DLA lie in the very noisy part 
of the spectrum affected by the Lyman limit break from the absorber at z = 2.92, and are 
thus unobservable. No radio and sub-mm observations exist either, so no information exists 
on the molecules from the DLA. We attempted to search for UV absorption features of CO 
AX 0-0, 1-0 , 2-0, 3-0, and 4-0 near 1544, 1510, 1478, 1447, 1419 A respectively, but no 
significant features were detected with a consistent set of velocity components. From the 
non-detection of the strongest of these features (CO AX 0-0), we estimate the 3 a upper 
limit to the CO column density as log Nco < 13.80, including the effect of both photon 
noise uncertainty and continuum determination uncertainty. This CO column density limit 
is comparable to the A^co detected in a QSO DLA by Srianand et al. (2008b), but much 
smaller than that in the log Nm = 22.7 DLA toward GRB 080607 (Prochaska et al. 2009). 

The non-detection of CO may seem surprising because of the high molecular mass 
fraction suggested by Fig. 9, and because a substantial H2 fraction may be expected at 
such a high Nm in the interstellar medium of the Milky Way. However, it is consistent with 
the low observed C I column density A^ci ~ 10^^ cm~^. This value of A^ci is close to the 
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regime where H2 becomes self-shielded, so A^Ha may range anywhere from ~ 10^ to ~ 10^^ 
cm~^ (see, for example. Fig. 22 and 27 in Welty & Hobbs 2001). Such a range in H2 column 
density is entirely consistent with the inferred upper limit on CO (log A^co < 13.8). We 
also note that the Krumholz et al. (2009a) scenario has some caveats (e.g., the assumption 
of 2-phase equilibrium and the exclusion of warm H I) which may not apply in this case 
(for example since there is clearly evidence of warm gas in this DLA). Such cases of low 
A^H2 at high A^hi are not common, but not unprecedented. For example, Welty, Xue, & 
Wong (2012) report sighthnes toward 9^ Ori C in the Orion Trapezium region and Sk9 in 
the southwestern part of the main SMC bar that also show low H2 at high A^hi, presumably 
because of strong radiation fields. Indeed, SMC sightlines with log A^hi ~ 22 show a wide 
range of the number fraction of H2 from ~ 10~^ to 10~^ (e.g., Tumlinson et al. 2002). Thus 
we conclude the molecular fraction is low in the DLA toward Ql 135-0010 because of the 
low metallicity, low dust content (see sec. 4.2, 4.3) and/or strong radiation field. It will be 
interesting to determine the molecular content of this absorber more accurately with higher 
S/N spectra in the future. 

4- 1-2. Dust Depletion vs. H I Content 

Watson et al. (2006) reported a DLA with log A^hi = 22.6 toward the afterglow of 
GRB 050401. They reported [Si/H] = -1.7, [Zn/H]=-1.0, [Fc/H] = -2.1, and [Cr/H] = 
-1.7. Compared to this, the DLA toward Q1135-0010 has about the same metallicity, but 
a smaller amount of depletion. In particular. Si is essentially undepleted relative to Zn 
([Si/Zn]= -0.04 ± 0.04), compared to [Si/Zn] =-0.7 in the GRB DLA of Watson et al. 
(2006). We note, however, that the measurements for the DLA toward GRB 050401 have 
large uncertainties ([Si/Zn] = — 0.7± 0.5), so the latter could be consistent with the [Si/Zn] 
for the DLA toward Q1135-0010 within the uncertainty. While [Fe/Zn] and [Cr/Zn] are 
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lower compared to [Si/Zn] in the DLA toward Q1135-0010, they are not as low as in the 
DLA toward GRB 050401. For the afterglows of GRB 050730 {zabs = 3.969), GRB 071031 
{zabs = 2.692), and GRB 080413A {zabs = 2.435), Ledoux et al. (2009) report log A^hi = 
22.10, 22.15, and 21.85, respectively. These latter absorbers show metallicities [X/H] of 
-2.18 ±0.11, -1.73 ±0.05, and -1.60 ±0.16 (where X = S, Zn, and Zn, respectively). The 
corresponding relative abundances are [Fe/X] = 0.06 ±0.06, —0.04 ±0.02, and —0.13 ±0.07, 
all of which are surprisingly low. Indeed, Ledoux et al. (2009) show that the majority of 
the GRB DLAs in their sample he at the high-Ani , low [X/Fe] end in the plot of A^hi 
vs. [X/Fe]. Fig. 10 shows [Fe/X] vs. H I column density (where X = Zn in most cases, S 
in a few) for QSO DLAs compiled in Meiring et al. (2009b), along with GRB DLAs from 
Prochaska et al. (2007) and Ledoux et al. (2009), along with the DLA toward Q1135-0010. 
It is clear that the DLA toward Q1135-0100 lies in the less-populated central right part of 
this plot (high H I column density, moderately large [Fe/X]). While the relative abundance 
[Fe/Zn] observed for this DLA is not the largest, it is in the upper half of the range of 
[Fe/Zn] observed for QSO DLAs at z > 1.5. 

4.2. QSO Reddening 

The SDSS magnitudes of Q1135-0010 are u, g, r,i,z^ 21.22 ± 0.08, 19.03 ± 0.01, 18.47± 
0.01, 18.25 ± 0.01, 18.25 ± 0.02, respectively. The A{g — i) color, relative to the median 
{g — i) for the SDSS composite QSO spectrum at Zem — 2.907 (adopted from Richards 
et al. 2001) is 0.48. If most of this apparent reddening arises in the DLA z — 2.2068, 
the implied E{B -V) = A{g - i)(l ± Zabs)'^''^ / 1-506 would be 0.08, assuming the SMC 
extinction curve. There are two other strong absorbers in this sightline at z = 2.696 and 
z — 2.924, but their Ly-alpha hues are much weaker (rest frame equivalent widths of 1.08 
and 1.13 A, implying log A^hi < 19). Thus these other absorbers are not even sub-DLAs), 
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and are not likely to contribute much compared to the DLA with log A^hi = 22.0. However, 
part of the inferred E{B — V) for the DLA toward Q1135-0010 arises from the fact that the 
strong DLA line falls within the g band and removes a large section of the QSO continuum, 
so the true reddening is likely even lower. 

The implied amount of low reddening in this DLA is comparable to that in most QSO 
absorbers at z ~ 2. For example, the reddening found by York et al. (2006) in a study 
of SDSS Mg II absorbers at 1 < z < 2 is < 0.09, while that found by for SDSS DLAs at 
2.2<z< 5.2 is -0.0017 ± 0.0022 (Frank & Peroux 2010) or < 0.03 (Khare et al. 2012). 

4.3. Extinction Curve Fitting 

We note that the SDSS spectrum shows no evidence of a 2175 A bump. Thus the 
dust in this absorber does not appear to resemble that in the Milky Way. To quantify the 
reddening, we fitted the SDSS flux-calibrated spectrum with the Milky Way, LMC, LMC 
super-shell (LMC2), and SMC bar extinction curves, following the procedure described in 
Kulkarni et al. (2011). The original spectrum of the QSO was taken to be the composite 
QSO spectrum adopted from Van den Berk et al. (2001). Next, for a given assumed 
extinction curve, the optical depth at any wavelength in the absorber rest frame was 
calculated as t(A) = 0.921 yl;^ = 0.921 yly^;^, where Ax is the extinction at wavelength A, 
Ay — Rv X E{B — V) is the extinction in the V band, and Ca is the extinction curve, all 
quantities being in the absorber rest frame. The expected spectrum was then calculated as 
a function oi E{B — V) and a flux normalization factor and matched to the dereddened 
observed spectrum with a least squares approach. The Fitzpatrick-Massa (FM) parameters 
for the average extinction curves and the Ry values for the Milky Way, Large Magellanic 
Cloud (LMC), Large Magellanic Cloud supershell (LMC2), and the Small Magellanic Cloud 
(SMC) bar from Misselt et al. (1999) and Gordon et al. (2003) were adopted. Emission 
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and absorption lines in the QSO spectrum were excluded while doing the fitting so as not 
to bias the continuum level. Fig. 11 shows the extinction curve fits to the SDSS spectrum 
dereddened for Galactic extinction, which is shown here box-car smoothed by a factor of 9 
for the purpose of display. The best-fit E{B - V) values for the MW, LMC, LMC2, and 
SMC bar extinction curves are -0.04, -0.03, 0.03, and 0.02, respectively, with reduced 
values of 6.07, 6.70, 6.20, and 5.30, respectively. Thus, among the four extinction curves, 
SMC bar appears to fit the best. In fact, all four extinction curves suggest low reddening, 
although none of these extinction curves fits well. (We note that the negative E{B — V) 
values for the MW and LMC extinction curves are not physically meaningful; they just 
indicate that there is no significant reddening, and in fact, there is no satisfactory fit for 
those extinction curves. Of course, as in previous studies utilizing similar techniques for 
extinction curve-fitting of QSO spectra (e.g., Srianand et al. 2008a), an individual QSO's 
spectrum could be different from the composite QSO spectrum template.) 

We thus conclude that the DLA toward Ql 135-0010 does not show much extinction. 
This is surprising compared to most interstellar clouds with high A^hi in the Milky Way, 
which show large extinction. The low extinction observed is consistent, however, with the 
lack of detectable CO hues, and low molecular fraction inferred above. 

4.4. Cooling Rate and Star Formation Rate Surface Density 

The absorber toward Ql 135-0100 offers a unique opportunity to estimate the star 
formation rate in the absorbing galaxy in two independent ways. As discussed earlier, the 
tentative detection of the Ly-a emission suggests a star formation rate of 9.6 ±2.1 M© yr~^ 
(although this detection needs to be confirmed with higher S/N spectra). An additional 
constraint on the SFR can be obtained from the detection of C II* A1336 absorption, 
following Wolfe et al. (2003). 
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Most of the interstellar cooling is due to the fine structure line emission of [C II] A 158 
fim, which arises from the ^^3/2 to transition in the ground 2p term of C II. The 
cooling rate can be expressed as Ic = Ncn-huui^ui/Nui, where Nqu* is the column density 
of the C II ions in the 2P3/2 state, while Aui and huui are the coefficient for spontaneous 
photon decay and energy of the ^^3/2 to ^Pi/2 transition (e.g., Pottasch et al. 1979). 
The UV transition of C II* A 1336 can be used to infer A'^cii*- For the absorber toward 
Q1135-0010, the range of A^'cii* consistent with the data is 1.01 x 10^^ < A^cii* < 4.66 x 10^^ 
cm~^. y This implies 2.6 x 10^^^ < h < 1-2 x 10^^^ erg s~^ per H atom. 

In Fig. 12, we plot the cooling rate vs. H I column density for interstellar clouds in the 
Milky Way, along with the corresponding measurements for DLAs from Wolfe et al. (2003) 
and the DLA studied here. The ISM measurements are adopted from Lehner et al. (2004) 
and shown separately for low, low+intermediate, intermediate, and high-velocity clouds 
in the Milky Way. Most DLA measurements seem consistent with the anti-correlation 
between 1^ and A^hi seen in the Milky Way interstellar clouds. The uncertainty in the C II* 
column density for the DLA toward Q1135-0010 mentioned above results in a fairly large 
uncertainty in the cooling rate. If the cooling rate for the DLA toward Q1135-0100 is close 
to the upper limit, it would be clearly very distinct from both the ISM clouds in the Milky 
Way and other DLAs. 

^This range of A^cii* corresponds to the C II column density range consistent with the 
C II A1334 line, i.e. 3.02 x 10^^ cm'^ < Ncu < 2.28 x 10^^ cm^^. Of course, we cannot 
rule out Ly-a forest blends, but it would be remarkable if the blends occurred exactly at the 
positions of the major velocity components expected in C II* from the other second ions. 
While we cannot completely rule out that some of the C II* absorption may be high-velocity 
C II, no other low ionization species in this DLA show absorption at such high velocities. 
Therefore this possibility seems remote. 
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Based on the observed gas- phase abundances [Fe/H] and [Si/H] (the latter almost 
equal to [Zn/H]), we estimate the dust-to-gas ratio to be about 0.076, using equation 7 
of Wolfe et al. (2003). Using this dust-to-gas ratio, the range of cooling rate Ic inferred 
from C II*, and comparing to the predicted coohng rate vs. gas density from calculations 
similar to Fig. 5 of Wolfe et al. (2003) for a range of SFR surface density ■0* values, one 
can constrain the cold neutral medium (CNM) density and hence ip^,. While a detailed 
calculation of the Ic vs gas density relation would give a more exact estimate of -0*, given 
the observational uncertainty in Ic, here we estimate the approximate range of ■0*. To do 
this, we make use of the fact that the observed Fe and Si abundances and lower limit on 1^ 
for the z = 2.2068 DLA toward Q1135-0010 are all within 0.1-0.2 dex of the corresponding 
quantities for the DLA at the comparable redshift z = 2.04 toward Q0458-02 from Wolfe et 
al. (2003). Thus, using Fig. 5b of Wolfe et al. (2003), we estimate that -2 < log V'* ^ 
Mq yr~^ kpc~^ from the limits on 1^ derived above. It is interesting that the SFR surface 
density expected from the observed H I column density on the basis of the Kennicutt (1998) 
relation V'* = [2.5 x lO""^] x [A^hi/1.26 x lO^^]^-^ would be in the middle of this range (log 
= -0.87). 

Adopting log ■0* ~ — 1 M© yr"-*^ kpc~^, and using the SFR 10 Mq yr"^ estimated 
in section 3.4 from the Ly-a emission flux along the line of sight, it may seem that the 
line-of-sight extent of the DLA absorber should be d ~ 10 kpc in the simple assumption 
of the star formation being distributed uniformly within the DLA. However, since both 
estimates of SFR are along the pencil beam sightline toward the QSO (with a projected 
beam size ~ pc), it appears that the star formation is highly clumpy (otherwise the implied 
DLA line-of-sight extent would be enormous). 
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4.5. Electron Density 

4.5.1. C II* Absorption 

The ratio of column densities of C II* and C II can be used to constrain the electron 
density in the gas, for an assumed gas temperature, because the C II upper level is expected 
to be populated by coUisional excitation, and depopulated by radiative de-excitation 
in the intervening DLA toward Q1135-0010. Equihbrium between these processes gives 
-/Veil* /-/Veil = neCi2(T)/^2i, where A21 = 2.29 x 10"^ s"^ (Nussbaumer & Storey 1981). 
The colhsion rate coefficient is given by Ci2(T) = [8.63 x 10~6Qi2/(c/iT0-5)] exp(-Ei2/A;T) 
(Wood & Linsky 1997), where gi — 2, £^12 = 1.31 x 10~^^ erg, and the colhsion strength ^12 
depends on temperature only weakly. At T~ 7000 K, Q12 = 2.81, giving C12 = 1.43 x 10~^. 
For the absorber toward Ql 135-0010, using the maximum C II and minimum C II* column 
densities (see section 4.4), we get Ncii*/-/Vcii > 4.4 x 10~^, which implies the electron 
density rie > 7.1 x 10~^ cm~^. A roughly similar electron density lower limit is obtained 
even if the gas is much cooler, say T — 500 K, since at this temperature Q12 = 1.82, giving 
C12 = 2.90 X 10^^, and hence rig > 3.5 x 10~^ cm~^. By contrast, if the minimum C II 
and maximum C II* column densities are used, we get the electron density rig < 24.7 cm~^ 
for T = 7000 K {n^ < 12.2 cm-^ for T = 500 K). Of course the upper limit is very high 
for either temperature considered, and the actual electron density is likely much lower (see 
section 4.5.2 below). 

We note that the maximum C II column density would imply a highly supersolar 
[C/H], reminiscent of the C-enhanced absorbers reported by Cooke et al. (2011) for much 
more metal-poor absorbers. The true C II (and C II*) column densities may be intermediate 
between the two limiting values discussed above, and the true [C/H] may not be as extreme. 
Correspondingly the electron density and cooling rates would be intermediate between the 
limiting values discussed above. 



4-5.2. Si ir Absorption 



Additional constraints on the electron density can be obtained from the Si II* 
absorption, which also appears to be present in some of the velocity components for 
this DLA. To our knowledge, this is the first detection of Si II* in an intervening QSO 
absorber, although Si II* absorption has been reported in GRB afterglows (e.g. Savaglio 
et al. 2012 and references therein). Fig. 13 shows velocity plots for the Si II* 1264.7, 
1265.0, 1309.8, and 1533.4 lines. (The Si II* 1265 panel shows the combined contribution 
of the A 1264.7 and A 1265.0 lines as the solid green curve and the contribution of Si II* A 
1265.0 alone as the dashed blue curve.) The green curves show our Voigt profile fits, with 
the main contributions coming from velocity components near -73, -31, and -19 km s"^. 
These fits yield a total column density (summed over all the velocity components) of log 
A^siii* = 13.70^°:°^, giving Nsiip/A^siiii = (1-62 ± 0.18) x 10"=^. The Si II*/Si II ratio can be 
used to estimate the electron density assuming equilibrium between coUisional excitation 
and spontaneous radiative de-excitation in this intervening DLA. Using the collisional 
excitation rate for Si II [given by = 3.32 x 10"^ (T/IO, 000)-°-^ exp(-413.4/r) cm^ s'^ 
e.g., Srianand & Petitjean (2000)] and the spontaneous radiative de-excitation rate for Si II* 
^421 = 2.13 X 10~^ s~^, the corresponding estimate of the electron density is Ug — 0.53 ± 0.06 
cm-3 for T = 500 K or 0.91 ± 0.10 cm'^ for T = 7000 K. 

4.6. Implications for Co-moving Density of Neutral Hydrogen 

The co-moving density of neutral hydrogen flni is sensitive to the presence of high 
A^Hi systems, since it constitutes an integral of the product of the H I column density and 
the column density distribution /(A^hi) which is roughly shaped as a Gamma function or 
a combination of power laws with negative slopes. We note, along the lines of Noterdaeme 
et al. (2009), that DLAs might have not been identified in the Sloan Survey because 
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the presence of the DLA trough significantly reduces the QSO flux and hence decreases 
the corresponding S/N of the spectrum. If many such high N(HI) column densities were 
missing, this would affect the shape of the column density distribution /(A^hi) and hence 
the resulting neutral gas mass, Qhi (Prochaska et al. 2005, Noterdaeme et al. 2009). Given 
the small redshift region (< 0.1) covered by the SDSS but not sampled by prior studies 
(e.g., Prochaska & Wolfe 2009; Ledoux et al. 2009), and the number density of strong 
DLAs, we expect very few if any other absorbers similar to the super-DLA among the 
SDSS DR7 data. Given the rarity of such systems, we anticipate that they will not have a 
significant impact on the estimates of ^hi- 

On the other hand, such high Ahi systems can provide unique insights into the nature 
of dust, metals, and molecules in the most gas-rich galaxies. It would be interesting to 
obtain radio / sub- mm observations of Ql 135-0100 for foUowup studies of molecular gas in 
this unique environment. It would also be of great interest to attempt to image the galaxy 
responsible for the super-DLA in this sighthne. We plan to pursue these directions in the 
near future. 

Note: Toward the very end of the refereeing process of this paper, we learned about 
another analysis of the same object by Noterdaeme et al. (2012), who used VLT X-shooter 
and Magellan MagE spectroscopy, together with our VLT UVES data. The Ahi value 
derived by Noterdaeme et al. (2012) (log Nhi = 22.10 ± 0.05) agrees closely with our value 
of 22.05 ± 0.10, as do the abundances of various elements (Si, Cr, Fe, Ni, Zn) that they 
derived from their X-shooter data and our UVES data. The dust reddening estimated by 
these authors from the X-shooter spectrum (0.04 for an SMC extinction curve) is also close 
to our own corresponding estimate (0.02). The Lyman-alpha emission is also confirmed in 
the X-shooter spectrum, and in fact shows a resolved double-peaked profile, giving a flux 
within a factor of 2 of our estimate from the SDSS spectrum (and within a factor of 1.25 of 
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the estimate from our VLT UVES spectrum) . 
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Fig. 1. — The SDSS spectrum of Q1135-0010 that enabled us to detect the super- 
DLA. The DLA trough is labeled. 
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Fig. 2. — The redshifted damped Lyman-alpha line from the VLT UVES spectrum 
of Ql 135-0010. Top: The solid black curve shows the un-normalized data, resampled 
to a dispersion of 0.42 A per pixel, (i.e., 30 times lower resolution than the original 
data), while the dashed black curve near the bottom indicates the 1 a uncertainty 
in this resampled spectrum. Bottom: The solid black curve shows the resampled 
spectrum after continuum normalization. The solid red and dashed blue curves 
overlaid on the data show the Voigt profile fits corresponding to log A^hi = 22.05 ± 
0.10, all centered at z = 2.2062. For comparison, we also show as dot-dashed orange 
curves the Voigt profile fits for log A^hi = 21.85 and 22.25, which are clearly too 
extreme. 
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Fig. 3. — Velocity plots of key metal lines in the Zabs = 2.2068 absorber toward Q1135-0010. The 
green curves indicate the multi-component Voigt profile fits. The dotted vertical lines indicate the 
positions of the velocity components. 
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Fig. 4. — Velocity plots of additional metal lines in the Zabs = 2.2068 absorber toward Q1135- 
0010. The green curves indicate the multi-component Voigt profile fits. The dotted vertical 
lines indicate the positions of the velocity components. The dotted and dashed vertical lines 
in the top panel show the positions of the velocity components for Ti II A A1910.612, 1910.954. 
In the bottom two panels, the dashed blue curves indicate the contributions from C II A 1334 
and O I A 1302, the remaining contributions coming from blends with C II* A 1336 or Ly-a 
forest lines. The C II A1334 panel shows the maximum C II contribution consistent with the 
data. The dot-dashed vertical lines in the C II A 1334 panel denote the contribution from 
C II* A1336. Unrelated lines of systems at other redshifts are blocked out in grey for display 
purposes. 
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Fig. 5. — Relative abundances vs. Zn II column density in different velocity components 
of the Zabs = 2.2068 absorber toward Q1135-0010. 
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Fig. 6. — Abundances of various elements relative to Zn, after summing over all the 
velocity components, plotted vs. condensation temperature. Three points are plotted for 
each element: the black circles denote the values, along with 1 a uncertainties, for the 
Zabs = 2.2068 DLA toward Q1135-0010, while the red squares and blue diamonds denote, 
respectively, the corresponding relative abundances for warm and cold diffuse interstellar 
clouds in the Milky Way. 
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Fig. 7. — Velocity plots of higher ionization species in the Zabs = 2.2068 absorber toward 
Q1135-0010. The green curves indicate the multi-component Voigt profile fits. The dotted 
vertical lines indicate the positions of the velocity components. For C IV, the dashed vertical 
lines indicate the positions of the velocity components in C IV A 1550. 
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Fig. 8. — The tentative Lyman-a emission feature seen near the center of the DLA trough, 
(a) Top: A closer look at the region near the bottom of the DLA profile. Shown are VLT 
UVES data resampled to 0.42 and 0.98 A resolution, with no further smoothing. Also 
overplotted is the unsmoothed SDSS spectrum, (b) Bottom: VLT UVES data resampled 
to 1 A resolution, with and without box-car smoothing. The vertical dashed lines indicate 
expected position of Lj-a emission at velocity separations of 0, +300, and -300 km s~^ 
with respect to z = 2.2068. 
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log N(HI) 

Fig. 9. — Metallicity (based on firm detections of Zn or S) vs. log A'^hi for QSO DLAs, GRB 
DLAs, and tlie DLA toward Ql 135-0010. The solid orange line shows the "obscuration 
threshold" suggested by Boisse et al. (1998). The remaining lines indicate the curves 
corresponding to molecular Hydrogen core "covering fractions" of 0.01, 0.05, 0.5, and 1.0, 
adopted from Krumholz et al. (2009a). The super-DLA toward Q1135-0010 is clearly 
well above the "obscuration threshold" , and consistent with H2 core "covering fraction" 
of ^ 0.5. 
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Fig. 10. — Depletion of Fe relative to Zn or S plotted vs. H I column density for QSO 
DLAs ai z < 1.5, QSO DLAs at z > 1.5, GRB DLAs at z > 1.5, and the DLA toward 
Q1135-0010. Error bars denoting 1 a uncertainties are shown where available. 
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Fig. 11. — SDSS spectrum of Q1135-0010 after dereddening for Galactic extinction, along 
with the best-fitting Milky Way, LMC, LMC supershell (LMC2), and SMC bar extinction 
curves. The SDSS spectra have been box-car smoothed here by a factor of 9 for display 
purposes. For each extinction curve, the intrinsic QSO spectrum, assumed to be the 
composite QSO spectrum from Vanden Berk et al. (2001), was reddened as a function of 
E{B — V) and flux normalization factor. The curves plotted correspond to the best-fitting 
parameters that give the least squares fit to the data. 
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Fig. 12. — Cooling rate inferred from C II* absorption line plotted vs. H I column density. 
The larger blue triangles show the estimated limits to the cooling rate we estimate for 
the DLA toward Ql 135-0010. The small black circles and triangles are detections and 
upper limits, respectively, for the sample of QSO DLAs in Wolfe et al. (2003). The small 
colored squares denote the measurements for low, intermediate, low -|- intermediate, and 
high- velocity interstellar H I clouds in the Milky Way compiled in Lehner et al. (2004). 
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Fig. 13. — Velocity plots of Si II* lines in the Zabs = 2.2068 absorber toward Q1135-0010. 
The green curves indicate the mult i- component Voigt profile fits. The dotted vertical lines 
indicate the positions of the velocity components. In the Si II* A 1265 panel, the dashed blue 
curve denotes the contribution of only Si II* A 1265.0, while the solid green curve denotes 
the total contribution of both Si II* A 1264.7 and A 1265.0 lines. 



Table 1: Voigt Profile Fitting Parameters for Various Velocity Components for Selected Ions in the z — 2.2068 DLA 



mean vel. 


b 


Cr II 


Zn II 


Fe II 


Ni II 


Si II 


(km s"^) 


(km s~^) 


(cm~^) 


(cm~^) 


(cm-2) 


(cm~^) 


(cm~^) 


-138 


9.4 


- 


- 


(1.35±0.17)el4 


~ 


(2.84±0.64)el4 


-115.5 


13.1 


(1.60±0.14)el3 


(3.07±0.86)el2 


(4.17±0.82)el4 


(5.31±0.64)el3 


(2.79±0.20)el5 


-75.5 


17.0 


(1.57±0.19)el3 


(4.73±1.23)el2 


(3.37±0.81)el4 


(6.36±0.87)el3 


(3.43±0.31)el5 


-60.0 


11.8 


(2.60±0.21)el3 


(4.57±1.00)el2 


(1.63±0.23)el5 


(5.86±0.84)el3 


(3.78±0.44)el5 


-39.0 


8.0 


(3.42±0.98)el2 


- 


(2.65±0.90)el4 


(1.43±0.44)el3 


(7.91±1.17)el4 


-16.5 


7.4 


(2.33±0.17)el3 


(1.33±0.14)el3 


(1.20±0.19)el5 


(8.65±0.90)el3 


(9.36±1.50)el5 


5.0 


8.3 


(1.22±0.12)el3 


(4.22±0.66)el2 


(3.48±0.96)el4 


(4.13±0.56)el3 


(2.16±0.21)el5 


30.5 


8.0 


(3.93±0.25)el3 


(1.08±0.10)el3 


(1.20±0.18)cl5 


(9.17±0.91)cl3 


(7.25±1.15)cl5 


62.5 


11.2 


(5.44±1.07)el2 


(1.07±0.54)el2 


(2.44±0.26)el4 


(4.55±0.58)el3 


(6.12±0.91)el4 


103.0 


12.4 






(1.58±0.43)el3 


(1.06±0.46)el3 


(5.50±0.47)el3 


145.0 


11.3 










(2.33±1.71)el2 


227.5 


15.5 






(2.24±0.53)el3 




(3.97±0.38)el3 



Table 2: Voigt Profile Fitting Parameters for Various Velocity Components for Additional Ions in the z = 2.2068 DLA 



mean vel. 


h 


C lit 


C IP 


oit 


Ti II 


P II 


Cu II 


(km s~^) 


(km s~^) 


(cm~^) 


(cm~^) 


(cm~^) 


(cm~^) 


(cm~^) 

V / 


(cm~^) 


-138 


9.4 






> 3.45el4 








-115.5 


13.1 


8.68el8 


5.57el4 


> 1.39cl5 


(8.3 ± 0.72)el2 


(1.74 ± 1.08)el3 




-75.5 


17.9 


7.56el8 


1.84el4 


> 8.82el4 


(5.28 ± 0.88)el2 


(3.48 ± 3.09)el3 


(4.97 ±3.37)el2 


-60.0 


11.8 


2.32el8 




> 2.58el5 




(1.22 ± 1.11) el4 


(5.22 ± 31.20)ell 


-39.0 


8.0 






> 6.55el4 


(9.69 ±4.85)ell 


(5.63±6.77)el3 


(1.33 ± 0.37)el3 




7 4 








(-\ c,4 _|_n c,«'ip-| 9 






5.0 


8.3 




1.04el4 


> 7.88el4 


(1.04 ±0.54)el2 




(1.65 ± 2.07)el2 


30.5 


8.0 


6.72el7 


1.05el4 


> 4.45el5 


(7.02 ±4.97)ell 


(3.57 ± 7.81)el2 


(4.88 ± 17.40)ell 


62.5 


11.2 




2.23el3 


> 5.88el3 


(1.22 ±0.80)el2 






103.0 


12.4 


5.00el4 


5.20el2 


> 2.59el2 






(6.84 ± 18.30)ell 


145.0 


11.3 


1.06el5 


9.58el2 








(6.81 ±17.90)ell 


227.5 


15.5 


5.30el6 


1.19el3 






(2.81 ± 1.40)el3 


(3.56 ±2.09)el2 



t The column densities for C II, C II*, and O I in some components could be higher because of hne saturation. 
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Table 3: Total Gas-phase Column Densities and Abundances in the z = 2.2068 DLA 



Element 


log (X/H)e 


log l^total 


[X/H] 


[X/Zn] 


C 


-3.61 


< 19.36 


< 0.92 


< 1.98 





-3.31 


> 16.30 


> -2.44 


> -1.38 


Si 


-4.46 


16.49 ±0.03 


-1.10 ±0.10 


-0.04 ±0.04 


p 


-6.54 


14.58 ±0.21 


-0.93 ±0.23 


0.13 ±0.21 


Ti 


-7.08 


13.28 ±0.04 


-1.69 ±0.11 


-0.63 ± 0.05 


Cr 


-6.35 


14.15 ±0.02 


-1.55 ±0.10 


-0.49 ±0.03 


Fe 


-4.53 


15.76 ±0.03 


-1.76 ±0.10 


-0.70 ±0.04 


Ni 


-5.78 


14.67 ±0.02 


-1.60 ±0.10 


-0.54 ± 0.03 


Cu 


-7.74 


13.56 ±0.10 


-0.75 ± 0.14 


±0.31 ±0.10 


Zn 


-7.37 


13.62 ±0.03 


-1.06 ±0.10 


0.0 
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Table 4: Voigt Profile Fitting Parameters for C IV, Si IV in the z = 2.2068 DLAf 



mean vel. 


^'SilV 


A^siiv 


^CIV 


iVciv 


(km s~^) 


(km s~^) 


(cm~^) 


(km s~^) 


(cm~^) 


-318 


4.6 


(6.16 ±3.39)el2 


4.6 


(6.55 ±1.53)el2 


-212 






19.6 


>1.51el3 


-202 


4.9 


>3.05el3 


4.9 


>5.55el3 


-176 


4.8 


>8.44el3 


4.8 


> 1.78el4 


-152 


13.3 


>1.07el4 


13.3 


>3.38el4 


-112 


13.1 


>3.51el3 


13.1 


>2.94el4 


-78 


17.0 


>9.79el3 


17.0 


>2.58el4 


-60 


11.8 


>2.93el2 


11.8 


>1.35el4 


-39 


8.0 


>1.26el4 


8.0 


>2.15el4 


-15 


7.4 


>2.07cl3 


7.4 


>2.35el4 


9 


8.3 


>1.18el4 


8.3 


>2.09el4 


34 


8.0 


>2.53el3 


8.0 


>1.49el4 


60 






17.7 


>5.12el3 


68 


11.2 


>2.03el3 


11.2 


>2.48el3 


105 


15.3 


>1.07el4 


12.4 


>3.26el4 


138 






6.3 


>1.48el3 


157 


11.3 


>9.76el2 


11.3 


>7.68cl3 


181 


13.3 


> 1.31el4 


13.3 


>7.96el3 


203 






7.6 


>1.32el3 



t The column densities for C IV and Si IV in most components could be higher because of 
line saturation. 



